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R. L. Jaffe* 

Department of Physics, Center for Theoretical Physics, and Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

Quantum mechanics predicts an exponentially small probability that a particle with energy greater 
than the height of a potential barrier will nevertheless reflect from the barrier in violation of classical 
expectations. This process can be regarded as tunneling in momentum space, leading to a simple 
derivation of the reflection probability. 



I. INTRODUCTION 

The most famous application of the semi-classical approximation in quantum mechanics is the calculation of the 

barrier penetration probability. P{E) = exp |— -| dxy / 2m(V(x) — ~E)\, for a particle of energy E to tunnel through 

a potential V(x), when the height of the potential exceeds the particle's energy for a < x < b. Elementary derivations 
of this result can be found in most textbooks [1]. Quantum mechanics also predicts that a particle with energy 
E, greater than a barrier V(x) for all x, will reflect from the barrier. This can also be estimated using semiclassical 
methods. A derivation can be found in Landau and Lifschitz's well-known presentation of the semiclassical mcthod[2], 
but the derivation is more complicated and less intuitive than the derivation of the barrier penetration factor. 

In this paper we show that reflection above the barrier can be understood as barrier penetration in momentum 
space, making the physical origin of the phenomenon more obvious and reducing the derivation to that of ordinary 
barrier penetration. This appears not to have been noticed in the past. 

An added benefit of this view of reflection above the barrier is that it leads to a simple derivation of the transition 
probability in the adiabatic approximation — the Landau-Zcner formula — again based only on the physics of barrier 
penetration. 

The paper is organized as follows: In Section II the problem of reflection above the barrier is shown to be barrier 
penetration in momentum space. In Section III the result is shown to be identical to the result derived in Ref. [2]. In 
Section IV a few simple examples are presented. In Section V the transition probability in the adiabatic approximation 
is recast as reflection above the barrier in time (as opposed to space) and the Landau-Zener result is derived. 

II. TUNNELING IN MOMENTUM SPACE 

To see that reflection above the barrier can be viewed as barrier penetration in momentum space, consider the 
classical motion of a particle incident from the left on barrier V(x), as viewed both in coordinate and momentum 
space. The familiar coordinate space depiction of the classical motion is shown in Fig. 1(a). For simplicity, take the 
barrier to be symmetric, V(—x) = V(x), so V max = V(0), and adjust the zero of energy so that V max = 0. The 
asymptotic value of the potential as x — > ±oo is — Vq. When E is negative, for example E < in the figure, the particle 
moves up the x-axis until it reaches the classical turning point at x = —xo(E < ), where V(—xq) = E < , then it turns 
around and retreats down the negative x-axis. Quantum tunneling gives a small probability that the particle makes 
a transition to x = +xq(E < ) and then propagates up the positive x-axis. When E is positive, for example E > in the 
figure, the particle propagates unimpeded up the x-axis from — oo to +oo. 

Now examine the problem from the point of view of momentum space, as shown in Fig. 1(b). For an energy E, 
the range of the momentum is bounded between — p max and +Pmax, where p m ax = \Jlm(E + Vo), as marked by the 
circles on the trajectories in the figure. For positive energy, for example the particle incident from the left in 
coordinate space moves down the positive p-axis to the minimum classical momentum, po(i?>) = \J2raE > , turns 
around, and retreats back up the positive p-axis. This corresponds to the classically expected motion over the barrier 
in coordinate space. Reflection above the coordinate-space barrier occurs because quantum tunneling gives rise to 
a small probability that the particle makes a transition to p = — po(-E>) and then propagates on down the negative 
p-axis. Having chosen Vmax = 0, the classically forbidden region in momentum space always ranges from — po(i?>) to 
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FIG. 1: Classical motion in the presence of a potential barrier in coordinate space (left panel) and momentum space (right 
panel). On the left the classical particle either reflects from the barrier when E = E < is negative, or passes over it when 
E = Ey is positive. On the right the classical particle reflects from the (always parabolic) barrier when E = i?> is positive 
and passes under it when E = E< is negative. 




FIG. 2: Phase space trajectories corresponding to barrier penetration (in blue) ABCD and reflection above the barrier (in red) 
EFGH y for the harmonic oscillator. 



Po(E > ) or vice versa. When E is negative, for example E < in the figure, the particle propagates unimpeded down 
the p axis from p ma x(^<) = \/2'm-(E < + Vq) to —p max (£'<), corresponding to reflection from the barrier in coordinate 
space. 

The analogy is particularly clear for the case of a barrier that is an inverse harmonic oscillator, V(x) — -~^ax 2 . 
In that case the classical orbits are hyperbolae in phase space as shown in Fig. 2. The barrier penetration trajectory 
(E < 0) is shown in blue. It progresses classically from A, at large negative x and positive p, to the classical turning 
point B, then jumps to C and continues on to large positive x and positive p. The refiection-above-the-barrier 
trajectory (E > 0) is shown in red. It progresses classically from E, at large positive p and negative x, to point F, 
jumps to G, and on to large negative p and negative x. The symmetry of the harmonic oscillator Hamiltonian under 
the exchange p <-> x (up to a change of scale) guarantees that both processes have analogous tunneling descriptions. 

To proceed quantitatively, consider the Schrodingcr equation for the momentum space wave function, (j)(p), 

(£ + y( ^ ) )^ )= ^ (P) ' (1) 
where V has been written as a function of the operator x — > ifi-r-, which is appropriate in momentum space. Un- 
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less V(x) is a sum of powers, V(ih-^) is an integral operator 1 . Things simplify significantly in the semiclassical 
approximation, where we expand the logarithm of the wave function in powers of H, 

4>(p)=vq>(la(j?) + <D(h )) . (2) 



H 

Substituting into eq. (1) and keeping the leading term in h, we find 



2- + V{x{jp))=E, (3) 



where x{p) = —<r'(p) = — ^2. This equation can be solved for cr(p) in terms of the inverse function V 1 , 

a{p) = J" 'dp'x(p') = - J" ' d V 'V-\E- V ^) (4) 



In classically forbidden region of momentum space, —\l2mE < p < y / 2mE, the inverse function, must be 

defined by analytic continuation from the domain £ < where it is defined. For example, in the case of the inverse 
harmonic oscillator, Vuo(x) = — \ax 2 , V^ {^) = ±\/— 2£/a, and 



Note that V~ x is always imaginary in the classically forbidden zone. The probability of reflection above the bar- 
rier, |![H(-E)| 2 , is the exponentially small probability of barrier penetration in momentum space, which we obtain by 
integrating J dpx(p) from — po(E) to po(E), and choosing the sign of x(jp) that corresponds to exponential suppression, 



which is our result. 



III. RELATION TO RESULT OF LANDAU AND LIFSCHITZ 

In their treatment of the semi-classical approximation, Landau and Lifschitz[3] derive a formula for the probability 
of reflection above the barrier that looks quite different from our eq. (6). They express the probability of reflection 
above the barrier as the exponential of a contour integral in the complex coordinate plane, 



|91ll(£)| 2 =cxp(--Im J dzi 



P(z) , (7) 



where p(z) = y / 2m(E — V(z)). zq is the complex value of z in the upper-half z-planc at which p(zq) = 0, and z\ can 
be any point on the real z-axis, since p(z) is real everywhere on the real z-axis (remember E > V(x) for all x). For 
a symmetric potential with maximum at x = 0, zq is on the imaginary z-axis, zq = iyo, so for convenience we take 
z\ = 0, and rewrite eq. (7) as 



/ 4 rva 
\9\ LL (E)\ 2 =cxpl [ --ImJ it 



dyp{iy) . (8) 



This result can be converted to ours using integration by parts, 



2/o rPo 

dy p(iy) = yp(iy)\o° + dpy(p). (9) 
o Jo 



1 If V(x) = !dqV{q)e-^/ h , then V(ih-^)<f>(p) = J dqV{q)e q ^ 0(p) = / dqV(q)<t>(p + q). 

2 It is possible to compute the 0(K) term as well, but as in coordinate space, it does not affect the exponential barrier penetration factor. 
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The surface term vanishes because p(iyo) = (This is how zq = iyo was chosen.). p(0) is the classical momentum at 
y = 0, so p(y = 0) = y/2mE = p , and the function y(p) is the solution to eq. (3) with x = iy: p 2 /2m + V{iy) = E, 
whence y = — iV~ 1 (E — p 2 /2m). Substituting into cq. (8) converts it to our result, eq. (6). So, in essence, the relation 
between the two results is summarized by J pdx = — J xdp. The derivation we have presented not only simplifies the 
mathematics, but also "explains" reflection above the barrier as barrier penetration in momentum space. 

IV. EXAMPLES 

A few simple examples may make the application of eq. (6) more transparent. We consider the inverse oscillator, a 
sech 2 x barrier, and a (1 + x 2 )^ 1 barrier. 

1. Vho(x) = -|ax 2 

The inverse of the harmonic oscillator potential is given by eq. (5), 

y H0 1 (£-i^ v / 2£-p 2 /™- (10) 
2m y/a 

Combining eqs. (6) and (10) we obtain a reflection probability of 



(2 pV2mE fly 
-T / dp\ -{E 



V 2 s \ 2irE 

f" ))=e~—, 11 
2m 



where uj = y/a/m is defined in analogy with the frequency of a normal oscillator. In this simple case, the exact result 
is known [5], 
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2wE 

e 



I^Ho(-E)|^ xa , 

1 + e~ — 

and agrees with the semiclassical estimate to leading order. 

2. V 2 (x) = V (sech 2 (x/a) - 1) 

In this case 



vy( E -pV2m)= ± .c,h- \J E _ p y° m + Vn )- M 

In the classically forbidden region, E > p 2 /2m, and the argument of the cosh -1 is less than unity, where cosh~ 1 (y) = 
icos -1 y. Substituting into eq. (6), we obtain the probability for reflection above the barrier, 



/ 2 f V2^E 

\X 2 (E)\ 2 =exp(-- dp 

\ n J-V2mE 



cos 1 



Vn 



y E -p 2 /2m + Vol I 
exp r_^V2^(v^T^-V^)) • (13) 



3. V 3 (x) = -V x 2 /(x 2 +a 2 ) 

In this case 



Vf\E-^/2 m) = iaJ v f + -^ 2m , iorp 2 <2mE, (14) 
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FIG. 3: (a) The function f(t) that parameterizes the time dependence in the adiabatic approximation, and (b) The time 
dependent eigenvalues of H(t). 



and the reflection coefficient can be expressed in terms of elliptic integrals, 

I*™* - sxp (d + ^) - 7K(^)) , (15) 

where 7 = Vo/E and E(x) and K(x) are complete elliptic integrals [6]. 

As E — > 0, the probability of reflection grows and the sensitivity to the details of the shape of V(x) disappears. 
This can be seen explicitly by taking E/Vq -C 1 in examples 2 and 3. In both cases the reflection probability reduces 
to the harmonic oscillator form with the parameter a replaced by the curvature at the top of the barrier, 7 — > 2Vo/a 2 
(the same for both the sech 2 and 1/(1 + x 2 ) barriers). 



V. THE TRANSITION AMPLITUDE IN THE ADIABATIC APPROXIMATION 



As a final remark, we use our approach to calculate the leading contribution to the transition amplitude in the 
adiabatic approximation. According to the adiabatic theorem if a Hamiltonian, H(t), changes slowly enough in time, 
a system initially prepared in an energy eigenstate will not make a transition to a different energy eigenstate, provided 
the two states remain well separated in energy. If the Hamiltonian changes slowly and smoothly, violations of the 
adiabatic theorem are exponentially small. For the important case of two levels that come close to crossing, the formula 
for the transition rate in the adiabatic approximation is due to Landau and Zcncr[7]. Qualitatively one can regard 
the Landau-Zencr transition as "tunneling" from one energy eigenstate to another. In light of the ideas developed 
in this paper, it is therefore natural to seek a description of the transition process in the adiabatic approximation as 
reflection above the barrier in the variable conjugate to energy, the time. 3 

For simplicity I consider two states and a Hamiltonian 4 

H(t) = (/W _^ } ) , (16) 

where f(t) is a function that decreases monotonically with t and is scaled by a parameter T, the "transition time", 
that determines the adabaticity. We assume that lim^ioo f(t) = =pE , with E 3> e, so near t = the levels come 
close to crossing, but in fact repel. A typical form for f(t) is shown in Fig. 3(a). 
At any time the eigenvalues of H(t) are 

E ± (t)=±Vf(t) 2 + e* = ±E(t), 



3 This idea is not new, indeed in Ref. [2] Landau and Lifschitz show that this problem is formally analogous to reflection above the barrier 
and obtain the transition probability by referring back to that problem. The argument presented here is quite similar but avoids the 
difficulties of analytic continuation to complex time. 

4 Without loss of generality we take f(t) to be antisymmetric, e to be real, and Tt H = 0. 
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and the eigenvectors are 



♦-«> - . (IT) 



where 



taaff(t) = 7 ^. 

A generic example of the energy level diagram is shown in Fig. 3(b). 

In the adiabatic limit the full, solutions to the time-dependent Schrodinger equation are simply proportional to 

&=(*), 8 

*±(t) =(j)±{t)e^S tdt ' E V\ (18) 

corresponding to E± = ±E(t). 

Suppose that at large negative time the system is prepared in the eigenstate 

#(*)-► ^ + (-oo)e-*^*/' l I as^->-oo, (19) 

where E = \JE 2 + e 2 . At large positive times the system will have a probability close to unity to remain in the state 
4>+ , but there will be a small probability to find the system in the state 0_ in violation of the adiabatic approximation, 

1(t>+(oo)e- lBt / h + m(t)-(oo)e +tBt/h , as t -> +oo . (20) 

|£H| 2 is the transition probability we seek to compute. Note that eqs. (19) and (20) define a problem analogous to 
scattering in one dimension. Identifying t with x, eqs. (19) and (20) describe a wave of amplitude T incident from the 
right (t 3> T), giving rise to a transmitted wave of amplitude unity moving off to the left (t -C — T) and a reflected 
wave of amplitude moving back to toward the right (t ^> T). 

To exploit this analogy we must convert the time dependent Schrodinger equation into a form that resembles the 
Schrodinger equation for scattering in one dimension. First write the time dependent Schrodinger equation for this 

problem without approximation. Let ^(t) = (^j^j ■ Then it$! = H^> implies 

iha = fa + eb 

ihb = ea - fb. (21) 

By differentiating the first equation and substituting the second, we obtain a second order differential equation for 
a(t), 

-h 2 'd- (f 2 + e 2 )a-ihfa = 0, (22) 

and an equation for b(t) in terms of a(t), 

b=-(iHa-fa). (23) 
e 

The equation for a(t) is a Schrodinger equation that describes scattering of a "particle" of mass 2m = 1 and energy 
e 2 from a "potential" 

V(t) = -f 2 (t)-ihf{t). (24) 
Eq. (23) merely fixes b(t) given a(t) and is irrelevant to our considerations. 



We do not need to allow for an adiabatic phase here because we are not considering a closed path in the space of states, and therefore 
the phase can be transformed away. 
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When T is large, f(t) is slowly varying, and the semiclassical method applies. The second term in V(t) (—ihf) 
is higher order in h and can be ignored 6 , and eq. (22) describes motion of a zero energy particle above the potential 
barrier V(t) — — (f 2 (t) + e 2 ). Without further ado, we can compute the transition rate by correspondence with the 
arguments leading up to eq. (6), 

V(x) - -f(t), 

E -> e 2 , (25) 
2m -> 1. (26) 



So the transition probability in the semiclassical approximation is 



l^l 2 = cxp (~ J" dp Tm{f-\y/j? 



e 2 )} , (27) 



which is our final result. 

Eq. (27) applies for an arbitrary slowly varying function f(t). The standard Landau-Zcncr formula applies to a 
linear level crossing, where f(t) = t/T. In this case / _1 (£) = T£ and 

= cxp (_HL f e dp ^T^\ = cxp (-H^.) , (28) 



which is the Landau-Zener formula. Of course eq. (27) applies to any slowly varying function of time. 
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This "quantum" contribution to the potential takes care of generating the prefactor cos 8/2 that modulates the semiclassical exponential 
in a(t). 



